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Abstract 



The kinetic Freeze Out process of a pion gas through a finite layer with time-like 
. normal is considered. The pion gas is described by a Boltzmann gas with elastic 

CN| \ collisions among the pions. Within this model, the impact of the in-medium pion 

mass modification on the Freeze Out process is studied. A marginal change of the 
Freeze Out variables temperature and flow velocity and an insignificant modification 
of the frozen out particle distribution function has been found. 

^ 1 Introduction 

One of the greatest discoveries in ultra-relativistic heavy-ion physics has been the creation of the 
Quark Gluon Plasma (QGP) at Conseil Europeen pour la Recherche Nucleaire (CERN) in 2000 
■ [TJ and at Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) 

in 2005 pj. In fact, there are compelling experimental signatures which transformed the QGP 
from a theoretical prediction into a precise observational science: elliptic flow V2, Jet quenching, 
Strangeness enhancement and Constituent Quark Number scaling can hardly be understood 
Q\ \ without the clear statement that this new state of matter has been achieved. However, during 

the last years it turned out that the produced new state of matter is more similar to a strongly 
coupled or strongly interacting Quark Gluon Plasma (sQGP), which has more characteristics of a 
liquid than of a weakly interacting plasma of quarks and gluons [3j HI [5j El El [91 [TO] ; for a recent 
comment on the term "sQGP" see [11J. While there is no doubt about a QGP phase transition 
of Quantum Chromodynamics (QCD) at T c = (173 ± 8)MeV [TO], it seems that a QGP of 
. freely moving partons can be reached only at higher energy densities and temperatures beyond 

the critical temperature T c . These new insights imply that further experimental signatures 
are certainly needed to understand not only the physical features of sQGP, but also how the 



new experimental facts do coincide with the predictions of the fundamental theory of strong 
interactions. More knowledge of this more complicated new state of matter are now expected 
from the intended experiments at the Large Hadron Collider (LHC) at CERN starting up very 
soon. In the following we will not distinguish between the terms QGP and sQGP, but want to 
keep in mind that the new state of matter is more complicated than expected from the early 
theoretical predictions. 

The evidence of a QGP cannot be proven directly. Instead, we have to trace from the 
observables at the detectors back to this very early stage of the heavy-ion collision. Obviously, 
the more accurate the description of the subsequent processes after forming the QGP is, the 
more accurate will be the picture and the understanding of this new state of matter. 

One promising theoretical method in this respect is the hydrodynamical approach based 
on the assumption of local thermal equilibrium. According to several theoretical studies, the 
produced QGP reaches a local thermal equilibrium very rapidly within (0.3 — 0.5) fm/c for gluons 
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and (0.5—1.0) fm/c for the quarks [131 13 IS] • Experimental data indicate a source size of less 
than 10 fm and less than 10 fm/c time extent. This strongly indicates a rapid pre-hadronization 
[T71 [18], which is also supported by the recent observation of constituent quark number scaling 
of collective flow data. Especially, when the expanding system reaches a temperature T < T c 
hadron states of high multiplicity, containing mostly pions, e.g. [19[ I20j. are formed. The pre- 
thermalization of quarks and quark clusters or pre-hadrons results in the local thermalization 
of pions |21l I22j . and even most of the low lying (s) hadronic states. 

Subsequently after or even during the hadronization, the chemical and thermal Freeze Out 
(FO) of the hadrons happens, where the hydrodynamical description breaks down and transport 
theoretical approaches are needed. First, the inelastic collisions among the hadrons cease, that 
is the so-called chemical FO at T c h- Immediately or simultaneously followed by the thermal FO 
at T t h, where also the elastic collisions among the hadrons are abondened. The FO process is 
essentially the last stage of the heavy-ion collision process and the main source for observables. 
An accurate FO description is therefore a basis for an accurate understanding of the initial states 
produced in ultra-relativistic heavy-ion collisions. 

A rigorous approach of the FO scenario from first principles is given by the Boltzmann 
Transport Equation which is a rather difficult assignment of a task. Even more, recently it 
has been recognized that the basic assumptions of Boltzmann Transport Equation are spoiled 
at the last stages of kinetic FO process |23[ I24j . and a more involved Modified Boltzmann 
Transport Equation has to be solved. The reason for that is because the characteristic lenght 
scale, describing the change of the distribution function, becomes smaller than the mean free 
path A at the last stages of the kinetic FO process. Thus, phenomenological models which can 
describe the kinetic FO process in a simplified manner by taking into account the main features 
of a typical FO process only, become rather important. 

Such a phenomenological description of kinetic FO process is usually modeled by two dif- 
ferent, in some sense even opposite, methods: A FO modeling through a hypersurface of zero 
thickness, and a FO modeling through an ininite space-time volume. Recently, the kinetic FO 
through a layer with finite thickness has been developed, for the case of space-like in [23] and 
for the case of time-like normal in |26j ; see also [27} [28] . This phenomenological approach makes 
a bridge between these two extreme FO models mentioned. So far, the impact of in-medium 
modifications of hadrons on the FO process has been considered only in Refs. [221 3Q[[3T]. In our 
investigation we will apply this recently developed FO model and consider a kinetic FO scenario 
through a finite time-like layer to study the impact of in-medium pion mass modification on the 
FO process. 

The paper is organized as follows: In Section [2] we give the needed basics of a transport 
theoretical description of a hot and dense pion gas. The kinetic FO process through a finite 
time-like layer is considered in Section The finite temperature mass modification of pions, 
embedded in a Boltzmann gas with elastic interactions among the pions, is examined in Section 
HI In Section[S]we present the results obtained, and in SectionOa summary is given. Throughout 
the paper we take c = h = ks = 1- 

2 Transport theoretical description of a pion gas 

Consider a system of N not necessarily conserved number of particles described by the one- 
particle distribution function f(x,p). This invariant scalar function is normalized by N = 
J d?r d 3 p f{x,p). Throughout the paper we consider a dilute pion gas where only elastic scat- 
terings among the pions are allowed, so that = (t,r) is the four-coordinate and p^ = (p°,p) 
is the four-momentum of the pion with p° = y/m 2 + p 2 . Then, the particle four-flow is defined 
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by, e.g. [33], 

d 3 p 
—o 

and the energy-momentum tensor is 



T» v = J ^■P U P V f(x,p). (2) 

The four-flow velocity of the medium can be defined as a time-like unit tangent vector at the 
wordline of the particles, i.e. = constant x (Eckart's definition). However, in case of 
non-conserved particles like the pions are, such a definition would not be convinient. Instead, 
for non-conserved charges or non-conserved particles the four-flow velocity of the medium is 
usually defined as a time-like unit vector parallel to the energy flow, i.e. ii M = constant x T^ u u u 
(Landau's definition): 

This tensor equation ([3]) is by definition valid in any frame. Obviously, in the rest frame of the 
gas, RFG, we would have u^ FG = (1,0, 0, 0) in Eckart's as well as in Landau's definition of the 
four-flow velocity. Any other frame of interest is related to RFG just by a Lorentz boost. 

From eq. (pQ) and eq. (|2|) one can define three linear independent Lorentz invariants: scalar 
particle density n, scalar energy density e and scalar pressure P, given by: 

n = N»u^ (4) 
e = UpT^Uv, (5) 

P = —T^A^. (6) 

where = g pu — u^Uu projects any four- vector into the plane orthogonal to the metric 
tensor g pv = diag(l, —1,-1,-1). We also note the invariant scalar entropy density, 

s = S»u„, (7) 

where the entropy four-current is defined by 

cPp 

Po 



S» = - [ ^-p»[f(x,p)lnf(x,p)-f(x,p)). (8) 
J Po 



In our investigation we will consider the FO process of an ultra-relativistic heavy-ion collision, 
and the hot region of the fireball shall be deemed to be in chemical (p n = 0) and local ther- 
mal equilibrium T{x). The conserved quantum numbers, {e.g. baryon number, electric charge, 
strangeness) are zero such that the thermal distribution can be characterized by the local tem- 
perature parameter T{x) of the fireball. Up to temperatures T < T c , most of the particles of 
such a system are the pions |19t 120] . which interact via elastic collisions with a cross section 
p.ciastic^ xhus, we will consider the invariant scalar functions eqs. - (J7J) for the case of an 
Boltzmann gas of pions where only elastic scatterings among the pions are allowed, and moving 
with a four-flow velocity u^. The particle distribution function in such a case is homogeneous 
f(x,p) = f(p) and given by the Jiittner distribution 

1 _____ P^ u v 



f cq {p,T) = g w ——exp[ ), (9) 
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where g n = 3 is the isospin degeneracy factor. Since the functions eqs. (jl]) - © are invariant 
scalars, they can be evaluated in any Lorentz frame. Especially, in the local rest frame RFG 
we obtain in case of an pion gas, characterized by the distribution function ([9]), the following 
expressions: 

n = ^mlTK 2 (a), (10) 

2 7T Z 

e = ^mlT [K! (a) + 3 K 3 (a)] , (11) 
P = J>^mlT*K 2 (a), (12) 

2 7T Z 

_ r 1 r? i 

(13) 



<?7r 2 
2 7T Z 



1 3 

TK 2 (a) + -77^^1(0) + -m w K 3 (a) 



where a = m n /T, and are the Bessel functions of second kind, see Appendix A. The Equation 
of State (EoS) P(n,T) of the pion gas follows from eq. (fTUj) and eq. (fT2|) . P = nT, and the 
thermodynamical relation Ts = e + P is also satisfied 0. In the limit of vanishing pion mass 
m n — > we obtain from eg. (Ilip and eg. (]12p the EoS of an ideal relativistic gas, e = 3 P. 

From these considerations we have seen that in RFG the particle density, energy density and 
pressure are only functions of temperature T. This implies that in any arbitrary Lorentz frame 
only two thermodynamical unknowns, temperature T and four-flow velocity u„, can enter the 
problem under consideration. To determine both unknowns, we need to have two differential 
equations, which can be deduced El from eq. ([3]) and eq. ([5]) [ 33] : 

da, = A *;; r / g , (15) 

de = u,dT^u v . (16) 
For the left side of eq. (1160 we obtain from the scalar invariant (|lip the following expression: 



de = ml [4 a K (a) + 8 K 1 (a) + 12 K 3 (a)]dT. (17) 

S7T Z 



These equations (|15p - (|17|) can be used to determine the temperature and four-flow velocity, 
while pressure P, particle density n and entropy density s would follow from their definitions 
in eqs. © and ([7]), respectively. The differential of the energy-momentum tensor needed in 
eq. (TL5j) and eq. (fT6j) follows from eq. ([2]), 

dT^ = J ^-^p u df(x,p), (18) 

according to which we still need a differential equation for the one-particle distribution function. 
This will be subject of the next Section. 



1 Recall, that a thermodynamical approach, by means of the canonical potential of an pion gas, 

d 3 p 

(2^)3 



n = -9- T I 7^ M 1 - cM-Po/T)) , (14) 



and with the aid of definitions of energy density e = —T 2 dil/dT, pressure P = — and entropy density 
s = dfl/dT, confirms the findings of eq. (fTTj) - eq. (fI3|) and the relation T s = e + P. 

2 For a proof of eq. (TT3)) see footnote 3 in [33], and for eq. (fib]) see also the remarks in the Appendix 
A of EH 
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3 Freeze Out Process within a finite time-like layer 



The scheme of an ultra-relativistic heavy-ion collision can be subdivided into three main stages 
characterized by their typical temperature parameter T: First, the initial stage at T c < T where 
a hot and dense parton gas is produced. Second, the stage at T prc „FO <T <T C where hadrons 
are formed. And third, the Freeze Out process at temperatures T post _FO < T < T pre _Fo where 
the hadrons freeze out. After the complete FO of the hadrons at T = T pos t_FO the particle 
interactions cease, i.e. the momentum distribution of the particles is frozen out and the hadrons 
move freely towards the detector. 

In this Section we are concerned with the third stage of the collision scheme, i.e. we start 
our investigation of FO process from the time of collision where the expanding system reaches a 
temperature T = T pre _Fo and the hadronization of the primary parton gas is considered to be 
completed. In the past, the FO process has been usually simulated in two extreme scenarios: a 
sudden FO on a hypersurface with zero proper thickness L = 0, or a gradual FO process during 
an infinite time and through an infinite space L — » oo. 

In this Section we present a model for a gradual FO through a finite layer, where the thickness 
L can be varied from zero to infinity, thus making a bridge between the two extreme schemes 
mentioned above. Within such a model the FO layer is bounded by two hyper-surfaces: the pre- 
FO hyper-surface with T = T pre _FO where the hydrodynamical description ends, and a post-FO 
hyper-surface with T = T post _FO where all the matter is frozen out. A covariant model of kinetic 
FO process within a finite layer has been recently developed, both for the case of space-like [25] 
and time-like [26] layer; see also [2"T1 |2"B]. 

In order to get an idea about the physical scales of the total FO time L of the FO layer, we 
recall that the collision time r co n between the pions, which are the dominant hadrons, depends 
on temperature r co n(T) = 12 f^/T 5 j35j [36l E7] , where /„- = 92.4 MeV is the pion decay constant 
(1 = 0.19733 GeV fm). At the pre-FO side of the layer there is a temperature of T prc _Fo ~ 
175 MeV and the collision time is small: r co n — l.lfm/c. The proper thickness L in time of the 
FO layer is taken typically of the order of a few (at least one) collision time at T = T pre _Fo ~ T c , 
i.e. L ~ (5...10)fm/c. 

Here, we will not repeat the theoretical developments of [251 GfH E7J [28] in detail, but should 
consider the essential steps relevant for our investigations. 

To describe the gradual FO process, the one-particle distribution function is decomposed 
into two components, an interacting part /, and a frozen-out part ff, 

f(x,p) = f l (x,p)+f f (x,p). (19) 

During the FO process the number of interacting particles decreases from pre-FO to post- 
FO side, where by definition the number of interacting particles tends to zero. As boundary 
conditions we assume on the pre-FO side of the layer a thermal equilibrium, i.e. a Jiittner 
distribution ([9]) for /, and ff = 0, while on the post-FO side fa vanishes; for an illustration see 
also fig. 1 in [23] . 

The space-time evolution of the interacting and non-interacting components during the FO 
should be modeled by the Boltzmann Transport Equation (BTE). We will use the relaxation 
time approximation and apply the escape rate P esc (x,p), describing the escape of particles 
from the interacting component fi to the non- interacting component ft. The FO is a strongly 
directed process, i.e. the gradient in one preferred FO direction, da^ = (dao,da), is much 
stronger than the changes in the perpendicular directions, we can neglect these FO gradients in 
the perpendicular directions. Then, the BTE can be transformed into the following differential 
equations [251 EE E7J EH] : 

<fcr M dpfifap) = -P esc (x,p) fi(x,p) + — [/eq(p) - fi(x,p)\ , (20) 
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da^ d tl f f (x,p) = P esc (x,p) fi(x,p) 



(21) 



Troughout the paper we shall use the notation 

d„f(x,p) = JLf( X)P ), (22) 

i.e. expressions like ()22p are not infinitesimal quantities. Note, the finite normal vector on 
hypersurface is normalized by da^ da^ = ±1 where upper sign is for time- like and lower sign 
is for space-like normal, respectively; note that da^ is also not an infinitesimal quantity but a 
finite vector (an explicit expression for time-like normal is given below). Here, p^ = (po,p) is 
the four-momentum of particle, and x^ = (t, r) is the four-coordinate of particle. The second 
term in eq. (|20p is the re-thermalization term (see below), which describes how fast the system 
relaxes into some thermalized distribution function / eq during a characteristic time scale rth- 
A Lorentz invariant expression for the escape rate is given by \24\ [25j [26"1 [2?] [25] . 



P^doji 

t L- x^ da^ p^ u p 



&P) = ~ r .».,_ 37^®(P^)> (23) 



where tq is the characteristic FO time. The — function is the Bugaev cut-off factor [35], which 
is important only for the FO in space-like direction. 

We have to insert the escape rate P esc into eqs. (|2(jp and (|21|) . In the following we will consider 
the FO process of a pion gas through a finite layer with a time- like normal da^da^ = +1, so 
that L becomes a thickness in time. We will work in the Rest Frame of the FO front, RFF, 
where da^ = (1,0,0,0). Thus, we obtain the following set of differential equations: 



r \L-tJ \p hl uv rth 



dtfi = ~- ( — ) ( -r^77 I fi + — t/eq(p) " M , (24) 



r \L-tJ \p^u^ 

Note again, the first term in eqs. (I24D and (I25p describes the transition of the pions from the 
interacting to the frozen out component. The second term in eq. (|24|) is the re-thermalization 
term [34|.I39| which describes how the interacting component relaxes to some thermal distribution 
/eq, where the parameters of it, T(t),u^(t), have to be calculated from the conservation laws. 
The strength of both terms are characterized by their typical time scales, the characteristic 
Freeze Out time To and relaxation time r t h, respectively [261 EH EH HQl ED E2J. 

In the case of fast re-thermalization r t h <C to, the interacting component can be choosen as 
equilibrated Jiittner distribution for all the times [25]. Then we obtain, with the aid of eq. (fT8|) . 
for the energy-momentum tensor of interacting component 

In the following we will give the components of eq. (|26|) in the RFF: 

dT?°(t,v,T,m n ) 1 L nT 1 



dt to L — t 4 7 u 

\ 0x (t,v,T, m n ) l dJf°(t,v,T,n 
dt v dt 



G^{m^v,T)-G+{m^v,T)\ , (27) 
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1 L nT b 2 



((3 + v 2 )K 2 (a)+aKi(a)) , (28) 



tq L — t 2 7 v 

dT xx (t,v,T, m n ) l dT° x (t,v,T, m n ) 

dt v dt 

T ( dNf{t,v,T,m^) 1 dN®(t, v, T, m„ 



+ - 



7« I dt v dt 

1 L nT 



a&fi(l + 3u 2 )if 2 (a) + &iTi(a)) , (29) 



tq L — t 2 

and the needed components of time derivative of the particle four-current are given by 
dN^(t,v,T,m n ) 1 L n 



dt Tq L - t 4 



-(G^(m^v,T)-G+(m^,v,T)) , (30) 



dN x {t,v,T,m^) ld/V^^m^) 1 L n ( AaK^a) 2a 2 K (ay 

Tq L — t 4 \ V V 



dt v dt 



(31) 



were we recall that a = m^/T and 6 = 70; the functions and K n are defined in the Appendix 
A. According to these expressions we also need the invariant scalar pion density n defined in 
eq. (jl]), and according to eq. (fl~5l) we also need the pressure P defined in eq. ([6]). Since they are 
invariant scalars, we can take the explicit expressions given in eq. (|lUp and eq. (|12p evaluated in 
the RFG for a Jiittner distribution. By inserting these results and eqs. (I27p - (|29D into eq. (I15p 
and eq. ()16p . we obtain a set of two differential equations for the two unknowns T and v. Taking 
""rff = 7(1)^0,0), u^ FF = 7(1, — v, 0, 0), duQ FF = ^vdv and dii^ FF = —j 3 dv, we obtain 
explicitly in RFF: 

£ = ^ (4ait ° (a) + SKM + 12A ' 3(a)rl 72 - + ° 2 T) • < 32 > 

dv 2vr 2 /1„ /x 3„,, L,,^ -1 / d7I po , 2 ^d7f E dT xx \ 

■ ^— r w — 7— • 

(33) 



, , -i^i(a) + -K 3 (a) + -K 2 (a) -v—*- + (1 + v , , . , 
dt c/ 7r m3TV4 w 4 JV 7 a v 7 1 dt eft dt 



Notice, that the degeneracy factor is actually cancelled against the same factor contained 
in the scalar particle density n of energy- momentum components, see eq. (jlUp . In the limit of 
vanishing pion mass — > the eqs. (|32H and (I33p simplify to 

dT 1 £_ 1 

dt" ~ ~Tq L-tl 7 ' 

(34) 

dv _ 1 L it; 

dt " r L-H 7 ' 1 j 

in agreement with the corresponding limit given in eq. (9) in [25] H The system of two differential 
equations ([32]) . ([33]) . together with the invariant scalars re, e, P in egs. (fTUj) - ([12]) and the 



3 Note, in case of a massless pion gas we have, due to Jiittner distribution, 

00 

47T 



.9tt 



(2tt) 3 




/" dpp 3 exp(-p/T) = 3^T 4 /7r 2 , (36) 
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components of energy-momentum tensor given in eqs. (|27p - f|29|) constitute a closed set of 
equations for the unknowns T(t) and v(t) inside the FO layer. Once these both unknows have 
been determined self consistently, all other quantities like particle density n(t), pressure P(t) 
and entropy density s(t) inside the finite layer can be determined by their expressions given in 
eqs. (HO]) - (fT3l). 

In our investigation, we are interested on the impact of in-medium pion mass modification on 
the FO process. Therefore, we have to implement a temperature dependent pion mass m n (T) 
in the given equations. The results of such an investigation can then be compared with the 
corresponding findings where a vacuum pion mass m^ or massless pions are implemented. We 
note, however, that in-medium modified pion mass not only changes slightly the distribution 
function, but also the pionic EoS and correspondingly the expansion dynamics, e.g. [43J. We 
note, that such a consideration could also be investigated, within the expanding model very 
recently developed in [97] . 

4 Pion Mass Shift at finite temperature 

Typical conditions inside the FO layer are high temperatures, typically between 100 MeV <T< 
175 MeV. Such extreme conditions imply a strong modifications of the hadrons in respect to 
their mass, coupling constants and decay rates are expected. One aim of our investigation is to 
evaluate how strong the impact of pion mass shift on kinetic FO process of a pion gas is. The 
kinetic FO process of a pion gas concerns elastic interactions among these particles and how 
these elastic scatterings cease. Therefore, we have to determine the mass modification of pions 
embedded in a pion gas with elastic scatterings among these particles. However, the impact of 
non-elastic interactions on mass modification of pions becomes relevant for a description of the 
chemical FO process. 

4.1 Pion mass in vacuum 

First, let us briefly re-consider the mass of a pion in vacuum, defined as pole mass of the pion 
propagator: 

II£(p) = i Jd A xe ipx (T w $ a (x)4> at (0)} 

= ipr^ > (38) 

p 2 — m w —Yj^(p) + ie 

where (0)0 = (0|O|0) is the vacuum expectation value of an operator O, a = 1, 2, 3 is the isospin 
index, T\y is the Wick time ordering, <& a is the second-quantized pion field operator and is 
the self energy of pion a in vacuum. The parameter m w is the so-called bare pion mass which 
enters the Lagrangian of the effective hadron model of QCD. The physical pion mass is defined 
as the pole of propagator eq. (|38|) . i.e. as the selfconsistent solution of 

ml = m n + ReZ n {p 2 = m 2 n ) . (39) 
while in [26] a Bose gas was assumed for the EoS: 

oo 

eBosc = .9tt 7^~)3 J d PP 3 i ex P(p/ T ) - 1] _1 = 9t, tt 2 T 4 /30 . (37) 
o 

The difference is only marginal. However, in eq. (9) of Ref. [26] we have actually to imply the relation 
^Bose = eBose/ (3 T) = g v tt 2 T 3 /90 valid for massless pions, to get the agreement stated above. 
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The vacuum pion mass is = 135.04 MeV, = 139.63 MeV. Despite enormous effort in 
Quantum field theory, a rigorous derivation of vacuum pion mass from first principles of QCD has 
not been found so far. That means the vacuum mass of any hadron can not been obtained from 
fundamental QCD without further assumption or new parameters. However, there are promising 
and sophisticated approaches which have provided some insights into this very involved issue. 
Among them there is Chiral Perturbation Theory, Lattice gauge theory, Current algebra, Dyson- 
Schwinger approach, Nambu-Jona-Lasinio model and QCD sum rules, which provide a link 
between the quark degrees of freedom of underlaying QCD and the hadronic degrees of freedom. 
Such approaches allow a derivation of the so-called Gell-Mann-Oakes-Renner (GOR) relation 

mm, 

mlfl = -2m q (M)o, (40) 

which allows to determine the pion mass in vacuum from the microscopic QCD quantities current 
quark mass m q and chiral quark condensate at x = 0, 

(M)o = ^ $u + <M>° • ( 41 ) 

Typical values are m q = (m u + m d )/2 = 5.5MeV, (fg) = -(245 MeV) 3 ; recall, f n = 92.4MeV 
is the pion decay constant. With these given numerical values, the GOR relation yields m n = 
138 MeV. The pion decay constant can be defined by 

<0|A-(0)|7r b (g)> = -iUq^Sab, (42) 

where the isospin indices a, b = 1, 2, 3 and 

A a x {x) = f(x) 7A476 ^*(x) (43) 

is the axial vector current operator, e.g. [67J; T a are the Pauli spin matrices, and iff = (u , d) T . 
From eq. (|42p we obtain the expression for the pion decay constant in vacuum, 

U = z^(0|i°(0)K(g)), (44) 

where, due to isospin symmetry, a can be chosen arbitrarily, e.g. a = 3 {i.e. there is no sum 
over index a). It should be noted that with the aid of LSZ reduction formalism [53\ 154] . the 
expression (|44p can be rewritten as vacuum expectation value over the scalar operator 

f n (0) = fd 4 yeM-iqy) (a y + ml)T w Aa(0)$^(y) , (45) 

Q J 

that means we basically have 

/* = (A)o- (46) 
4.2 Pion mass at finite temperature 

Now let us consider the case of pions in medium. At finite temperature the pions move in a hot 
and dense bath of hadrons. Accordingly, the in- medium pion propagator reads 



II£(p,T) = % Jd A xe ipx (T w $ a (x)6 a t(0)) T 



2 1 • (47) 



p 2 - m n -E«(p,T) + ie 
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The thermal Gibbs average (O)t of an operator O is defined by 



( )t = — —, — -, ~„ = —55 ; > ( 48 ) 

Tr (exp(-/3H)) £ ( n |exp(-/3ff)|n) 

n=0 

where /3 = 1/T. In general, for temperatures T >T C quark and gluon degrees of freedom would 
have to be included, i.e. the sum in eq. (1351) runs in general over quark and gluon eigenstates and 
over hadron eigenstates of H. However, in our investigation the temperatures are T < T c , i.e. the 
operator H is here an effective Hamiltonian describing the hadron system under consideration. 
Accordingly, we consider |n) as hadron eigenstates of H, and the sum runs over the spectrum of 
all these hadron eigenstates \n). Furthermore, throughout our considerations we are interested 
on Gibbs average of operators at x = 0, i.e. we always have (C')t = {0{x = 0))t- As in vacuum, 
the pion pole mass at finite temperature is defined as selfconsistent solution of 

ml (T) = ml + Re E w (p 2 = m 2 (T) , T) . (49) 

From this equation it becomes obvious why there is an in-medium modification of pion mass: just 
because the self energy ^ n (p, T) is now temperature dependent which changes the position of pole 
mass compared to vacuum. Such a change of the pole position is caused by elastic interactions 
and non-elastic interactions among the particles. To find an expression for the temperature 
dependence of the pion mass, we will follow a similar way as in vacuum, i.e. we will apply the 
in-medium GOR relation. Especialy, it is a well-known fact, that the GOR-relation continues 
to be valid also at finite temperatures, at least up to order 0(T & ) [351 1471 1351 1391 [501 ISTl 152] . 

m 2 (T)/ 2 (T) = -2m q $q) T . (50) 

We will take relation (|5(jp as a given fact, which would allow a determination of the temperature 
dependence of pion mass; the current quark mass m q as a fundamental parameter of QCD 
is of course independent of temperature, while the constituent quark mass M q is temperature 
dependent. For that we would have to know the temperature dependence of the chiral condensate 
and of the pion decay constant. That means, in generalization of the vacuum expectation values 
in eqs. (I4ip and (I46p . we have now to determine the expressions 

J2 (n\qqex.p(-/3H)\n) 

(W)T = 2^; , (51) 

E (n\exp(-pH)\n) 

n=0 

U(T) = ^ ; , (52) 

E (n\exp(-PH)\n) 

n=0 

where the operator f n is defined in eq. ([351) . Both of these expressions, (qq)T and f n (T) = (f^T, 
have been evaluated by means of several approaches. The aim is to find an expression both 
for chiral condensate and pion decay constant at finite temperature and consistently within 
the kinetic model description. Especially, in our approach there are only pions. Accordingly, 
the Gibbs average of an operator O runs over the diagonal pion states only, i.e. (7r|0|7r), 
(7r7r|0|7T7r), ... [611 162j. but does not include the diagonal matrix elements of heavier hadrons 
like Kaons, (K^lOlK^), etc. Such an approximation can be justified, since at low temperatures 
T <T C , pion states dominate the thermal average, while states containing heavier mesons with 
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a mass m n > m n are weighted with their corresponding Boltzmann factor ~ exp(— m n /T), i.e. 
they are exponentially suppressed [121 [2D]- For instance, in Ref. [121 EH] contributions of heavier 
mesons to the chiral condensate at finite temperature have been studied, where only marginal 
corrections were found: 5 percent corrections at T = 100 MeV and 10 percent corrections at T = 
150 MeV. The Gibbs average (118]) can be further approximated by one-pion states [58 ] [59 1 160]. 

(6) r = (O) + E / (2^2 po fr»l<V(P)> ex P ('Pi + ° ( T4 ) ■ (53) 



n=l 



The contributions of next higher order 0(T 4 ) are considered in Appendix C, according to which 
we will neglect multi-pion states in the temperature region considered, see also Refs. |61|, 162] . 
The pion states are normalized by 



7T 



! (pi)k m (p 2 )} = 2p° 1 (27Tfd^( Pl -p 2 )5 nm , (54) 



with isospin indices n,m = 1,2,3. With the aid of eq. ([53]) and according to eqs. ([5T]) and ([5"2~]) . 
we obtain 

^ P /_n/„\i4-^i_n/,^\ I P I ,/-, /-,-n /gg\ 



(Wh = mo + E / (27r / 2p o <*"(p)lf$l**(p)> exp + O (T 

U T ) = U + E / (2^0 <->)|i> a (s)vr>)> exp (-^j + O (t 4 ) 



(56) 



where in the last line the LSZ reduction of pion state n a (q) has been transformed back just 
after the thermal Gibbs average. The pion matrix elements can be determined using soft-pion 
theorem [58] ESI EH EH ES] by means of which we obtain (see eqs. ([71]) - ([77jl in Appendix B) 

(Mh = (M)o (i - ~ b 1 ( ! y))+o (r 4 ) , (57) 



UT) = /.(l-^^if^) ]+0fT 4 ) , (58) 



1^! r flZh 
12/| 1 v 

1 T 2 _ /m 



™*C0 = m^l + A^^^jj+o^ j (59) 

for the next higher order 0(T 4 ) see Appendix C. In order to derive relation (I59p . we have 
inserted egs. ([57]) and ([58]) into GOR relation at finite temperature eq. ([50]) . as well as the GOR 
in vacuum eq. (jlU)) has been applied. The function is @ 

oo 

g /■ g g 

BAz) = — ^ / dx V x 2 — z 2 exp(-x) = — =■ z KAz) , Yiva. BAz) = — ^ . (61) 



4 In case we would have taken a Bose distribution in eq. (|53[) . then in egs. ([57)) - ([5^)) the function £?i 
would have to be replaced by the function B 2 : 



oo 

B 2 (z) = % I dx yfx 2 - z 2 -\ , lira B 2 (z) = 1 . (60) 

■k 1 J exp(a;J — 1 z— o 

z 

However, the difference for pion mass shift using B\ or B 2 is marginal in the temperature region consid- 
ered, see fig. [U 
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Figure 1: The temperature dependent pion mass m n (T): Thick solid linepoints show the 
pion mass for a Boltzmann gas according to eq. (|59[) . The straight dotted line represents a 
constant pion mass in vacuum, = 138 MeV, according to the GOR relation in vacuum 
( 140]) . The dahed line represents the pion mass shift in case of a Bose distribution, i.e. 
when in eq. (159|) function B\ is replaced by function B 2 . The difference between a pion 
mass shift of a Boltzmann gas and of a Bose gas are marginal in the temperature region 
considered. The dotted line represents the result of ChPT according to eq. allowing 
a comparison of our results with ChPT. 



The results f)57[) - ()59[) are consistently valid for a pion gas approximated by a Boltzmann gas, 
i.e. with elastic interactions among the pions. Here, we will require the applicability of eqs. (j57H 
- (|59p in a temperature region T < T c , where the distance between pions is still large lj, and 
where the correlation between pions is still small [711 172"] . 

According to eq. (I59p the pion mass increases with increasing temperature, in line with other 
theoretical investigations. Especially, more sophisticated approaches yield very similar or even 
same results for the pion mass shift, e.g. Chiral Perturbation Theory (ChPT) [73], Nambu- 
Jona-Lasinio-model [73], QCD sum rules [75l [76], Linear Sigma Model [77J EH [721 EQ] , Mean 
Field Approximation [81] and Virial Expansion [82] , For instance, we can compare our findings 
(|57p - (j59[) with the results of ChPT, since at order 0(T 2 ) the pions are treated as free particles 
within ChPT, e.g. [56]. For (f<?) T [13 [55], for f„(T) p>] and for m w (T) [Ml [56], the ChPT 
yields following expressions: 




5 For instance, at T — 150 MeV the mean free path of pions is A ~ l/(n er°^ stlc ) ~ 2fm. 
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m^T) = 




The difference between our results given in eqs. (|57|) - (|59|) and the results of ChPT given in 
eqs. (j62j) - (j64"|) simply consists in the function B±, reflecting the fact that according to the above 
mentioned references the ChPT considers a Bose gas and the limit Bi — > 1. The result given 
in eq. (|62ft has also been derived in [57] and confirmed later on within the Sigma- model [68] . 
The given result in eq. ([63]) is confirmed, e.g. in [6TJ [68l [691 EQ] ■ In this respect we should refer 
the interested reader to Ref. [7U] where a Pade approximation of f n (T) has been established, 
valid for arbitrary temperatures T < T c . Nonetheless, to apply that result f n (T) one needs 
a Pade approximation for (jqq)T as well in order to be consistent within the entire framework. 
We also refer to [5H [52] where a stronger decrease of the pion decay constant f n (T) and a 
stronger increase of pion mass m n (T) for high temperatures were found within the framework 
of Dyson- Schwinger Equations. 

It is worth mentioning, that for in-medium pion mass shift there is an experimental result 
available |83[ [84"[ [85[ 186]. namely for the case of a deeply bound pionic state with a ir~ at finite 
baryonic density. Especially, a significant pion mass upshift of 19.3 ± 2.7 MeV in 207 Pb, i.e. 
approximatly at nuclear saturation density no — 0.17fm -3 , has been found. 

Another important experimental fact concerns the Constituent Quark Number scaling re- 
cently found [871 [88] . Such a scaling refers to an evident dependence of hadron elliptic flow V2 on 
the number of constituent quarks in the hadron under consideration. This experimental fact can 
be understood by assuming that FO and pre-hadronization after QGP occurs at an intermediate 
constituent quark stage [891 E0] • Accordingly, the pre-hadrons are made of constituent quarks 
which have a mass of about M q ~ (200 — 300) MeV. For the pre-pions (this argumentation is not 
valid for all the other hadrons, because they are considerably heavier than the pions) in QGP 
it would imply a mass much higher than their vacuum mass m n ~ 138 MeV. From this point of 
view a mass increase with increasing temperature and baryon density is in agreement with these 
experimental facts. Indeed, QCD sum rules [39] and QCD Dyson- Schwinger Equation [511 [52] 
predict a very strong pion mass up-shift near T c . In this respect we note, that a quark clustering 
like qq,qq,gg,qg etc. in the sQGP state has also been proposed in [91~1 19"2~]. 

5 Results and Discussion 

In this Section we will present and discuss the solution of the coupled set of differential equations 
(|32p and f|33|) for massive pions, once with the vacuum pion mass and once with a temperature 
dependent pion mass according to eq. (I50|) . In order to ascertain the impact of pion mass itself 
on the FO process, we will compare these results with the case of massless pions. 

The differential equations (|32|) and (|33|) are solved with the aid of Runge-Kutta-Method 
[93]. The limit — > given in eqs. (I34p and (I35p has been solved independently using MAPLE 
[94] . i.e. used as independent verification. The two boundary conditions of these first-order 
differential equations are the initial temperature on pre-FO hypersurface, T pre _Fo — T c = 
175 MeV, and the initial flow velocity on pre-FO hypersurface of the finite layer, u pre _Fo = 0.5 c. 
Both of these are typical values; note that such initial temperatures at small baryonic densities 
can be reached for instance at SPS [95]. In the numerical evaluation we have taken tq = lfm/c 
and L = 10 tq. The results are plotted in figs. M and El 

Let us first consider the temperature of the interacting component plotted in fig. [2j The 
temperature decreases in time, because the number of interacting particle decreases in time, see 
also fig. [H The difference between the case of a temperature dependent pion mass m n (T) (solid 
line) and the case of massive pions with constant vacuum mass m n (dashed line) is marginal. The 
reason for that can be understand by means of the small left figure of fig.® At high temperatures 
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Figure 2: Temperature of interacting component as function of time, when the matter 
crosses a finite Freeze Out layer. In the numerical evaluation we have taken tq = lfm/c 
and L = 10 r . Three situations of Freeze Out are considered: a pion gas with a tem- 
perature dependent pion mass m 7r (T) (solid line, top), the case of pions with a constant 
vacuum pion mass m n (dashed line, middle), and for massless pions m n = (dotted line, 
bottom). As the system expands, the temperature of interacting component drecreases in 
all three cases, since the particles with larger momentum freeze out faster. The lighter the 
pions are, the faster the Freeze Out proceeds. In the left small figure the marginal impact 
of pion mass shift even at extremly high temperatures is shown. The right small figure 
shows the fast drecrease of temperature at the post Freeze Out surface. The substantial 
difference between the curves for massive pions and massless pions shows how important 
the impact of pion mass on the Freeze Out process is. 



T ~ T c the pion suffers a mass up-shift of about 5 percent of it's vacuum mass, see fig. [TJ And 
indeed, at such high temperatures we find a difference of about 5 percent in the temperature 
slope between the case m n (T) and m n . However, at such high temperatures the impact of a 
pion mass is small anyway, because the averaged momentum of pions due to thermal motion is 
considerably higher than the pion mass itself. At lower temperatures the impact of a pion mass 
on FO process is stronger because the averaged momentum of pions becomes comparable with 
the pion mass itself. However, already at moderate high temperatures of about T ~ 135 MeV the 
pion mass up-shift is less than 5 percent, see fig. [TJ and we cannot expect any longer a significant 
impact of the mass shift on the FO process. In fact, while the in-medium modification of pion 
mass shift has almost no impact, it turns out that there is a significant deviation between the 
case of massive pions and massless pions. This is due to the fact mentioned, that at moderate 
high temperatures the pion mass becomes very comparable to the averaged momentum of the 
pion. In the right small figure implemented in fig. [2] we can see the fast dropping of temperature 
at the end of the FO process. This is because of the factor L/(L — t) in front of eq. (j32j) which 
becomes infinite at the post-FO side of the hypersurface. 
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Figure 3: Velocity of interacting component as function of time, when the pions cross 
a finite Freeze Out layer. In the numerical evaluation we have taken tq = lfm/c and 
L = 10 tq. Three cases are considered: massive pions with a temperature dependent 
pion mass m n (T) (solid line,top), massive pions with a constant vacuum pion mass m n 
(dashed line,middle), and massless pions m n = (dotted line, bottom). The line of in- 
medium modified pion mass (solid line) is above the line of a constant pion mass. This is 
due to the fact that the pion mass is a bit larger in the former case. During the Freeze Out 
process the fastest particles leave the interacting matter component first, so the velocity 
of the remaining non-interacting component decreases. The left small figure shows the 
tiny impact of pion mass shift, while the right small figure shows the fast dropping of 
velocity at the very end of the Freeze Out process. The remarkable difference between 
the curves for massive pions and massless pions shows the importance of pion mass on 
the Freeze Out process. 



Let us now consider the flow velocity of the interacting component plotted in fig. [31 The flow 
velocity decreases in time, because the number of interacting particles also decreases, see fig. [U 
As in case of temperature, the difference between a FO scenario with a temperature dependent 
pion mass m^T) (solid line) and a constant vacuum pion mass m n (dashed line) is negligible. 
The reason for such a behaviour is the very same as in case of temperature: At pre-FO side 
of the layer, there is a remarkable pion mass up-shift of about 5 percent compared to vacuum 
pion mass, which implies roughly a 10 percent impact on the flow velocity, as seen in the small 
left figure of fig. El However, at such high temperatures T ~ T c the pion mass is considerably 
smaller than the averaged pion momentum due to thermal motion, so that the impact of a pion 
mass is negligible. As the system cools down, the pion mass approaches rapidly the vacuum 
pion mass, see fig. [2j That means, at moderate temperatures of T ~ 135 MeV where a pion 
mass starts to have some impact of the FO process, the pion mass up-shift is already less than 
5 percent, so that the impact of in-medium modification of pion mass on FO process becomes 
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negligible. However, there is a considerable difference between the case of massive pions and 
massless pions, as can be seen in fig.\5\ The right small figure implemented in fig. shows that at 
the post-FO side of hypersurface the velocity falls down rapidly because of the factor L/(L — t) 
which becomes infinite at times near L. 

The fig. H] shows the rapid decrease of particle density of the interacting component. More 
than 90 percent of the particles get frozen out before a time scale of t ~ 3 To is reached, both 
for massive and massless case. 



0.25 




Figure 4: Density of interacting component as function of time, when the pions cross a 
finite Freeze Out layer. In the numerical evaluation we have taken To = 1 fm/ c and L = 
10 To- Three cases are considered: massive pions with a temperature dependent pion mass 
m^-(T) (solid line), massive pions with a constant vacuum pion mass m, (dashed line), and 
massless pions = (dotted line). During the Freeze Out process the particles leave 
the interacting matter component, so the density of the remaining interacting particles 
decreases. For the strong decrease of pion density we recall that for low temperatures 
the particles decreases exponentially: lim<r-^o n (T) = g n / '(47r 2 )m^T y / 27r/aexp( 
m^/T. 



-a): a 



In heavy-ion collision experiments the relevant experimental parameters are of course not the 
temperature T(t) or flow-velocity v(t) of interacting component. Instead, via the experimentally 
accessible spectrum of the particles dN/d?p one measures the post-FO distribution, which is 
formed in the applied model at the outer hypersurface, /fo(Px) = ff(t = L,p x ). The Cooper- 
Frye formula [96] 



Po 



dN 
d 3 p 



dVpo fpoip) ~ Po /fo(p) 



(65) 
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Figure 5: Left pannel: Massless pions. Solid line: FO distribution function /fo(Px) = 
ff(t = L,p x ), evaluated by means of eq. ( |25l) with T(t) and v (£) evaluated previously, see 
figs. (j2J) and (J3]). Dashed line: A fit of the above solid line by a thermal (i.e. Jiittner) 
distribution with T = 125 MeV and v = 0.49 c. Both distribution functions peak sharply. 
At low momenta p x , there is a considerable difference between a thermal distribution and 
the FO distribution function, i.e. the FO distribution function fpo{p x ) evaluated within 
the kinetic model applied is obviously not a thermal distribution. This is actually one 
of the main results of the kinetic FO model. Right pannel: Massive pions. Solid line: 
FO distribution function fpoiPx) = ff(t = L,p x ) evaluated by means of eq. ( 1251) with 
T(t) and v{t) evaluated previously, see figs. (jSJ) and (j3J), and with an in-medium pion 
mass m n (T). Dashed line: A fit of the above solid line by a thermal Jiittner distribution 
with T = 140 MeV and v = 0.4c and m n (T). There is a considerable difference between 
a thermal distribution and the FO distribution function for a massive pion gas at low 
momenta p x , a result which is in agreement with the corresponding statement in the 
massless case. The distributions for massless and massive pions peak at different momenta 
p x \ This is expected to effect flow measurements. The sharp peak is smeared out after 
summing up for fluid elements with different pre-FO flow velocities. 



allows us to calculate the particle spectrum (da^ = (dV, 0, 0, 0) is the time-like normal vector 
of an infinitesimal element of the post-FO hypersurface). In the one-dimensional model used, 
the volume V of the pionic fireball is ill-defined 0, so that only a qualitative comparison with 
experimental data, for instance the transverse spectrum, is possible, i.e. j^- ~ /fo(Px), where 
p x is locally pointing out in transverse direction. Due to the marginal impact of in-medium pion 
mass modification it is meaninful to compare the case of massive pions m n (T) with the case of 
massless pions, and both of each with a simple thermal distribution. 

In fig. [5] the frozen out distribution functions fpoipx) for massless and massive pions are 
plotted and compared with a thermal distribution. It is shown, that at low momenta the 
distribution function fpo(Px) evaluated within the kinetic FO model differs considerably from 
a simple thermal distribution, both for massless and massive pions. This confirms a very recent 
FO evaluation where a gradual FO with Bjorken expansion and with a constant pion mass 
m n = 0.138 MeV has been considered [HTJEH]- Notice, that the peak position of distribution 
function for the case of massive pions differs from the peak position for a massless pion gas. 
This is expected to effect experimental flow measurements. 

6 Since the volume of the pionic fireball is not defined, the global energy and momentum conservation 
cannot be applied. However, global conservation laws can lead to different results if and only if the 
curvature of the surface is large, thus leading to a significant divergence. Such a situation has been 
considered in |32j . where also related issues concerning an expanding system are discussed. 
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Figure 6: Contur plot of FO distribution function, fpoiPxiPy) = ff(t = L,p x ,p y ), evalu- 
ated by means of eq. (1231) with T(t) and v(t) evaluated previously, see figs. ([2D and ©• 
Solid line for the case of massive pions m^(T), Dashed Line for the case of massless pions. 
There is no curve for the case of massive pions if /fo — 1CF 2 because the FO distribution 
function in case of massive pions is smaller than this value, see fig. [5j It illustrates the 
importance of pion mass for an accurate FO description. The non-thermal asymmetry of 
the post-FO particles is a strong and dominant feature, even for a time-like normal FO 
hypersurface or layer. 



In fig. [6] the contur plot fFo(Px,Py) for the massive case m n (T) and the massless case m, = 
is compared. For large momenta p there is only a slight difference between the FO of massive 
and massless pions. However, at low p there is a significant deviation between a massive pion 
gas and the massless case. The figs. [5] and \5\ demonstrate that the mass of pions is relevant for 
an accurate FO description. 

6 Summary 

The strong interacting Quark Gluon Plasma is not directly accessible via experiments. In fact, 
one has to trace back from the experimental data to the initially formed new state of matter 
during an ultra-relativistic heavy-ion collision. This implies a detailed and accurate description 
of all subsequent stages of heavy- ion collision process. Within our FO model of a pion gas 
described by a Boltzmann gas, we have considered the problem how strong the impact of an 
in-medium modification of the pion mass on the kinetic Freeze Out process is. 

We have considered the kinetic FO process through a finite layer for a pion gas with finite 
pion mass and with massless pions. Throughout the investigation, all calculations have been 
consistently performed within the model of a pion gas with elastic scatterings among them and 
approximated by a Boltzmann gas. 

First, we have found, that there is a strong impact of the finite pion mass on the FO process 
compared to a FO process with massless pions. This result is highlighted in figs. [2] and [3l where 
a significant modification of the basic thermodynamical functions temperature T(t) and flow 
velocity v(t) of elastic interacting component of pion gas inside a finite FO layer has been found. 

Hereafter, we have investigated how strong the impact of in-medium modification of pion 
mass on the FO process is. To determine the pion mass at finite temperature, the generalized 
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GOR relation at finite temperature has been applied as a given fact, and the needed expressions 
(qq)t, fn(T) have been evaluated by a thermal average over one-pion states. The needed pion 
matrix elements were evaluated by the soft-pion theorem. Especially, the modification of effective 
pion mass is increasing with increasing temperature, becoming significant around 100 MeV and 
reaches ~ 5 MeV at T ~ 180 MeV, see fig. [D 

It has turned out, that the impact of in- medium pion mass shift m n (T) on the FO process 
is marginal compared to a massive pion gas with a constant pion mass m n . The physical reason 
for that can be understood as follows: A significant pion mass shift occurs at high temperatures 
T ~ T c . However, at such high temperatures the impact of the pion mass itself on FO process is 
negligible because the momentum of thermal motion is sizebly larger than the pion mass. During 
the FO process the temperature decreases rapidly, and after the initial ~ 2 To time it was below 
T ~ 100 MeV. But already at moderate temperatures T < 135 MeV, where the impact of pion 
mass on FO process becomes relevant, there is no longer a remarkable in-medium modification 
of pion mass. 

For a qualitative comparison with experimantal data and in order to investigate the impor- 
tance of a finite pion mass for the FO process, in figs. [5] the FO distribution function /fo(Px) 
has been plotted, both for massive m n (T) and massless = pions. The FO distribution 
functions have been compared with a thermal distribution. It is shown that at low momenta 
there is a considerable difference to a thermal spectrum, both for a massive and massless pion 
gas. The distributions peak at different momenta p x which is expected to effect experimental 
flow measurements. The contur plot in fig. [6] elucidates a remarkable difference between massive 
and massless pions at low momenta. Both figs. [5] and [6] show the importance of a finite pion 
mass for an accurate description of the FO process. 

The numerical order of in-medium pion mass shift at finite temperatures used in our model 
calculation is in some sense also encouraged from experiments which have found a pion mass 
up-shift of 19.3 ± 2.7MeV around baryon saturation density n ~ 0.17fm -3 [Ml El IBS]- 
However, we have to take into account, that the in-medium pion mass shift might be much more 
pronounced at high temperatures T ~ T c than used in our model calculations. For instance, 
QCD sum rules [19] and QCD Dyson-Schwinger Equation |51l [52] predict a strong pion mass 
up-shift near T c . 

Such a strong in-medium modification of pion mass would be also in line with the argu- 
mentation of quark pre-clustering inside the QGP for temperatures higher but near the QCD 
phase transition at T > T c , because the quarks will have a constituent quark mass of about 
M q ~ 200 — 300 MeV , so that the mass of pre-pions for such temperature regions would be much 
higher than the vacuum pion mass. Our results indicate that an earlier pre-hadronization and 
FO at higher temperatures would increase the sensitivity on the changing of effective pion mass. 

From experimental side, as mentioned in Section d such a conclusion is mainly supported 
by the experimental fact of Constituent Quark Number scaling, i.e. the dependence of hadron 
elliptic flow 1)2 from the number of constituent quarks of the hadrons. Therefore, for a more 
comprehensive analysis how strong the impact of in-medium modification of the pion mass for 
an accurate description of a heavy-ion collision process is, further insides of (pre-)pion mass shift 
near (and beyond) the critical temperature T c , both from theoretical as well as experimental 
side, are mandatory. Especially, the pion mass modification caused by non-elastic scatterings 
among the pions have to be evaluated in order to describe the chemical FO process. 

In summary, in our model study we come to the conclusion that, while a pion mass has 
a significant impact on the kinetic FO process, the temperature dependence of the pion mass 
m n (T) is negligible for an accurate description of the kinetic FO process. 
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Appendix A 

The function G^ (n = 1, 2) is defined by 

G±(m n ,v,T) = ^2 / dp p {j p 2 + m 2 ^ Ei ^ p 2 + ml 

o 

where Ei is a special case of incomplete Gamma-function [99J 



"fvp 



(66) 



Ei (a;) = / dt t~ l e~* . 



(67) 



In the massless case we get the temperature independent limits 

lim Gf = 2 1 1 



lim Gt 



3 7 3 {l±vf ' 
3 1 1 



~ 2 7 4 (1 ± v) A ' 

The function K n is the Bessel function of second kind [99], defined by 

2 n n! 



(68) 



K n {z) 



(2n) 



z~ n dx e~ x (x 2 - z 2 )"- 1 / 2 . 



Appendix B 



In this Appendix we will evaluate the matrix elements in eqs. (|55jl and (15611 with the aid of soft- 
pion theorem [58\ [63j [Ml [65"1 ]66j (a derivation of soft-pion theorem for the here used simplier 
case of "free" pions can be found in Appendix of Ref. jlUUj ) given by 



lim lim{7T m (p 2 )\d(x)\7T n ( Pl )) = -1<0| 

P2^° PI— *0 Jn 



2>n f 



0(x) , Q\ 



|0>, 



(70) 



where Q\ is the (time independent) axial charge, i.e. the spatial integral over zeroth component 
of axial vector current (|43|>: Q\ = J d 3 r Aq(t, t), and = Qa> an d n,m = 1,2,3 are the 
isospin indices. [A, J3]_ = AB — BA is the commutator. Here, we are interested in pion matrix 
elements of the following two-quark operators at x = 0: 



(71) 



Oo 



i=l a,j8=l 



3 4 



75^) = r X) X! (7o7/,75)a/3 (u'J - d'J d 1 ^ , (72) 



i=l a,/3=l 
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where the greek letters a, (3 denote Dirac indices, and i is the color index of quarkfields. The 
equal-time anti-commutators of the full QCD quark fields read for quark operators of the same 
flavor 

&(ri,t) , $(r 2 ,t)]+ = <5( 3 )(n - r 2 ) 5 a/3 & , (73) 

where [A, 5]+ = AB + BA is the anti-commutator, while quark fields of different flavor anti- 
commute. Using identities like [A,BC]- = [A,B]+C — B[A, C]+ and the anti-commutator 
relation (|73jl we obtain 

lim hm(7r m (p 2 )|fg|7r" ( Pl )) = * < |f$|0) <J™ , (74) 

J-K 



p 2 ->0 Pi— >0 



1 {3 -I 



lim lim (7r i (p 2 )|^|7r' 1 ( g )7r 1 (p 1 )) = - (0\A^(q)) = i - Qfl , (75) 

lim lim <7r 2 (pa)|^Jl7r 3 (g)7r 2 (pi)) = --^(0|ij k 3 (g)) = i g„ , (76) 

lim lim <7r 3 (p2)|4j|7r 3 (g)7r 3 (pi)) = 0, (77) 

p 2 ^0 Pi->0 



where in the last step of eqs. (|75j) and ([76]) we have used relation (|42 



Appendix C 



In our model we have taken into account the Gibbs average over one-pion states of an operator 
C, see eg. (I53p . Here we will consider the next order, i.e. the Gibbs average over two-pion states 
of an operator O: 



\ I n^n / 77735 n„n ex P I -777 ) exp 

n,m=l 



°( T4 ) 2! ^ / (2vr/L? / 



(2tt) 3 2p? J (2vr) 3 2p° 



x<7r m (p 2 )7r>i)|6|7r>i)7r m (p 2 )> , (78) 

where the factor (2!) _1 in front of (|78p circumvents a double counting by integrating over the 
permutations p\ «-> p 2 . With the aid of soft-pion theorem we obtain the pion matrix elements 
for chiral condensate and pion decay constant: 

]T lim lim(vr m (p 2 )7r"(p 1 )!fg|vr"(p 1 )7r m (p 2 )} = -|j <0|fg|0> , (79) 

n,m=l p 2->0 P1 ^0 Jtt 

£ lim lim(7r-(p 2 )7T> 1 )|i 3 |7r 3 (g)7r"(p 1 )7r m (p 2 )) = ± (0|i 3 |vr 3 (g)) = -i^ ^ , 

(80) 

where in the last line we have applied eg. (|42p . By means of these expressions we obtain 

<f$>T = ^)o(i-^5 1 + I ^^^+o(r 6 ), (si) 

m,(T) = m„ (l + 1| 2* + ^| 5?) + OGI*) , (83) 
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and B\ = B^m^/T). A numerical evaluation shows that the terms of order 0(T 4 ) contribute 
at most a few percent compared to the terms of order 0(T 2 ) in the temperature region we 
are interested. Within ChPT approach very similar statements are found, but it should be 
mentioned that in ChPT non-elastic particle interactions yield additional contributions to order 
0(T 4 ), while eqs. (|8T|) - (i83j) are the results for a pion gas with elastic particle collisions only. 
The dynamical non-elastic interactions among the pions change not only the given numerical 
values in eqs. (I8ip - (I83D but also the sign of the coefficients in front of the given order 0(T ), 
e.g. Ref. [62], where the result (|8ip is also mentioned. But we note again, that in ChPT 
the given coefficients of the order 0(T 2 ) remain untouched even when taking into account the 
contributions of non-elastic pion scatterings. 

Furthermore, in a pion gas with elastic interactions and to all orders in temperature we 
obtain by iteration 



that means all higher orders are factorial suppressed even at temperatures near T c . Again we 
underline, that eqs. (183]) - (f86|) are useful to study the impact of pion mass shift on kinetic FO of 
a pion gas (cease of elastic interactions among the pions). However, they are not applicable for 
modelling the pion mass shift impact on chemical FO process of a pion gas (cease of non-elastic 
interactions among the pions), because then the non-elastic interactions among the pions will 
change significantly the given series of higher order. 




(84) 




(85) 




(86) 
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